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Abstract 

. We show that a variety of complex processes can be viewed from the unified standpoint of 

^ ^ . scale-free nonlinear conservative cascades. Examples include certain turbulence models, perco- 

C I lation, cluster coagulation (aggregation) and fragmentation, 'coarse-grained' forest fire model of 

^ ■ self-organized criticality, and scale-free network growth. We classify such cascades by the values 

of three indices, and show how power-law steady spectra may arise. The power-law exponent is 
proven to depend only on the values of the three indices by a simple algebraic formula. 
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The cascade idea is central to our understanding of an amazingly broad scope of natural 

S'lenomena. One of the earliest mathematical cascade models, the Galton- Watson process 
, was suggested in 1874 as a solution to the "problem of the extinction of families." This 
kind of branching cascades was later used for theoretical treatment of chemical and nuclear 
chain reactions [2] , cosmic radiation , and statistics of polymers among other physical 
and biological problems. 

Cluster fragmentation is another common type of cascades [^. Typical examples of 
: ragmentation are rock fracture, breakup of liquid droplets and polymer degradation 
3]. Fragmentation preserves the total mass of clusters. Both branching processes and 
fragmentation are linear cascades, as the agents are assumed to not interfere with each 
other. 

Studies of turbulence have provided a wealth of examples of nonlinear cascades. Among 
them are the Richardson cascade of energy in the inertial range that underlies the now 
classical approach of Kolmogorov 0, the Kraichnan reverse enstrophy cascade for two- 
dimensional turbulence 0], and various random cascade models sug gested as theoretical 
explanations of the unsolved problem of turbulence intermittency |lO]. 

More recently the cascade idea has been frequently evoked in connection with a vari- 
ety of complex processes. A self -org anized branching cascade process was put forth as a 



.ea„-fle,d theory fo. aval^ch. U- A cascade .odel was suggested as a We-g.ained- 

approximation to the forest fire model, displaying self-organized criticality The classic 
result for the asteroidal size distribution [l3| has been re-considered from the perspective 

n n 

of collisional fragmentation cascades 114^. Mesoscopic rainfall n5 , cascades of reconnecting 

n n 

magnetic loops in solar flares [16||, and cascades of data networks [17( are further examples 
from this vast list. 

We remark, however, that other important complex systems and processes - e.g. multi- 
plicative ra„do,„ processes, percolafon Q. c.ust^ coagulafon M. Scheidegger mode, of 
river networks, directed Abelian sandpile model |20|], and scale-free random networks j2J| - 
can be advantageously considered as nonlinear cascades. 

Nonlinear cascades are characterized by strong interaction between the agents. Unlike 
linear cascades, they are are poorly understood. Their studies have usually relied on system- 
specific features and assumptions and, as a result, the general picture has remained unknown. 

The purpose of the present paper is to study a general class of scale-free nonlinear conser- 
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TABLE I: Some scale-free conservative cascades with stationary power-law spectra, n(s) oc s 
discovered in various contexts. For all of them t = 1 + (m + a)/h. 
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Directed Abelian sandpile model 
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Takayasu (1989) [25j 


Asteroid collisional fragmentation 
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Scale-free collisional fragmentation 
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vative cascades, and to show that seemingly different systems reveal common features when 
considered from this unified viewpoint. In particular, we derive a nonUnear cascade equa- 
tion, and demonstrate how it can give rise to stationary power-law spectra. We show that 
the power-law exponent depends only on three indices characterizing a particular cascade: 
the 'conservation law' index m, the scale homogeneity index a, and the nonlinearity index 
h. This result is illustrated by Table 1. 

Let us set the stage. We will think of cascades as statistical collective systems of in- 
teracting agents of different type. Examples of agents are turbulent eddies, coalescing rain 
droplets, avalanches, connected network clusters and polymer chains, genetic family trees, 
etc. As a result of an interaction of two or more agents, 'new' agents form, whereas 'old' 
agents annihilate. For instance, the result of hydrodynamic interactions between turbulent 
eddies leading to exchange of their energies, is viewed as annihilation of 'old' eddies and 
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creation of 'new' eddies having altered energies in their place. Any such interaction can be 
written formally as 

A,, + A,, + ... ^ A4 + A4 + ... ; or ^ u^As. = 0. (1) 

Si 

Here s labels various agent types, and can be discrete or continuous. In the latter case 
summation has to be substituted by integration. Here we will assume s to be a real number, 
and refer to it as 'size'. i/<j-'s are taken positive for 'new' agents and negative for 'old' ones. 
Although this notation is routinely used for chemical reactions, it is quite unusual in other 
fields, e.g. turbulence. 

We assume that agents of the same type are indistinguishable in statistical sense, and 
that a broad distribution of agent types emerges as a result of interactions. 

A cascade state is specified by the agent distribution, also called spectrum, n(t, s). Com- 
mon examples of agent distributions are the turbulent energy spectrum E{t^ k) in homoge- 
neous isotropic turbulence, where the scalar wavenumber k plays the role of s; the avalanche 
size distribution n(t,s); and the distribution of masses of clusters formed due to cluster 
coagulation or/and fragmentation. 

The primary goal of a cascade theory is to predict the spectrum n. We derive a novel 
formal nonlinear cascade equation, which describes temporal evolution of n, from the master 
equation [27[: 

Pt[N{s)] = J2 w[N{s),N^{s)]P^[N,{s)]-w[N,{s),N{s)]Pt[N{s)]. (2) 

Ni{s) 

Here Pt[N{s)] is the time-dependent probability of the microstate N, square brackets em- 
phasize functional dependence; w[N, Ni] is the probability of transition from the microstate 
A^i to the microstate per unit time; summation is carried over all microstates compatible 
with the interactions (Q). Then the macroscopic spectrum n{t,s) can be obtained as an 
average over all microstates: 

n(t,s) = (Ar(s))^5^Ar(s)P,[iV]. (3) 

N 

After multiplying (j21) by N{s), summing it over all microstates A^, regrouping the terms, 
and introducing AA^(s) = N{s) — Ni{s), we obtain 



(t, s) = ^AN{s)J2HN + ^N,N]Pt[N]. 



AN N 



When an interaction is of the type (^J), AN{s) has the form: 

AN{s) = u{s, {s,}) = ^sAs - s.). (4) 

Si 

Clearly, i/(s, {sj}) is the number of agents of size s created due to an 'elementary' interaction 
of agents of sizes {si}. Now we can define the interaction probability per unit time: 

W{{s,}, [Nis)]) = w[N{s) + u{s, {s,}), N{s)], (5) 

and its average: 

T{t,{s,})^{W{{s,},[N])). (6) 

The physical sense of T is interaction strength between the agents of sizes {st}. Thus, we 
arrive at the nonlinear cascade equation: 

n{t,s) = {{u{s,{s,}))), (7) 

where the following notation has been introduced: for any function /({sj}). 

The summation in (jH)) is carried over all sets {si}. At this point equation ((Tj) is not closed, 
as the interaction strength T depends on the microdistribution A^. The necessary closure 
can be provided for particular types of interactions by mean field-type approximations or 
renormalization procedures among other methods ^28] . However, we will show that for con- 
servative scale-free cascades important conclusions about asymptotic spectra can be found 
by-passing the closure problem. 

Let us illustrate the above ideas with two examples. For cluster coagulation jlflj ]: 

W{t; Si, Sj, Si+j] [N]) = K{si, Sj)Ns^Ns., for any i,j. 

Here K{i,j) is the probability that two cluster of sizes Sj and Sj aggregate per unit time. K 
is commonly called coagulation kernel. In the mean-field approximation: 

T = K{si, Sj)nsinsj, 
5 



and ((Tj) becomes the celebrated Smoluchowski coagulation equation: 

^ S— 1 OD 
Si = l Si = l 

As another example consider two-dimensional isotropic turbulence. In the inviscid limit 
the energy balance equation for a turbulent flow in a cyclic box of side D reads 

1 r°° 

dtE{k) = T{k), T{k) = -j^ dpdqTik,p,q), 

where 

T{k,p,q) = 27rA;Im{(D/27r)^(27r/|sin(p,q)|) 

X {kmSij + kj6im) {u*(k)Uj{p)Ura{q))}, 

(k = p + q, /c = |k|, p = IpI, g = |q|). 

This problem can be easily written in the form (|7j): 

Afc + z/pAp + z/gA, = 0, 

iy{s;k,p,q) = 5{s - k) + iyp5{s - p) + iyq5{s - q), 

Vp = T{p,k,q)/T{k,p,q), Ug = T{q,p,k)/T{k,p,q), 

T{k,p,q) = T{k,p,q). 

The condition that a cascade is scale-free has several consequences: first, the reaction 
type defined by (0) is independent of the physical scale of agent sizes. More exactly, any 
two sets of agent sizes, {sj} and {s-}, such that s- = Asj for some positive A, allow reactions 
of the same type, i.e. Ug. = Vgi . 

Next, interactions should have no characteristic size. In mathematical terms: if = 
Asi, N' = X-^N, then 

{wm,m)) = ^''{w{{s.},[N])). (11) 

For instance, for homogeneous isotropic turbulence a = 1 and this value can be traced to 
the gradient operator in the Navier-Stokes equations. For coagulation, a is the homogeneity 
degree of the coagulation kernel. 

Finally, the interactions must be scalable in agent density. Thus, if s- = Sj, A^' = AA^, 

{Wm,[N'])) = X'{W{{s,},[N])). (12) 



For example, for homogeneous isotropic turbulence h = 3/2, because the interactions depend 
on the third order velocity correlation function, and velocity scales as square root of energy. 
For coagulation, h = 2 due to binary collisions. 
Let the conservation law read: 

i^{sj, {si})sj' = 0, for any set {sj. (13) 

For instance, the energy conservation law for turbulence yields m = 0, while enstrophy 
conservation for 2D turbulence results in m = 2. Mass conservation in coagulation and 
fragmentation cascades yields m = 1. 

The nonlinear cascade equation can be used for treating open systems, e.g. a constantly 
stirred turbulent flow, a coagulating system with a source of the smallest clusters, or a driven 
self-organizing system. Let us assume (without limitation of generality) a direct cascade, 
i.e. s > So, where sq is the smallest size in the system. Then, we supplement equation (jZj) 
with the 'boundary condition': 

n(t. So) = no for all t, (14) 

and look for stationary solutions. We are particularly interested in the limit s ^ so- For 
scale-free conservative cascades this problem can be easily solved in a general way. Because 
the conserved quantity flux, E, is well-defined in the limit s/sq ^ oo, one can expect that 
solutions of (|7j), (fT^ depend on E and 'forget' the microscale so- These expectations are 
corroborated below. 

It follows from (|7j) and the definition of E that: 

E = Y.( W{{si}, [N]) ) ^s^sfe^S - Si), for any so<S<oo. (15) 

{Si} Si 

Here 9{x) is the step function. A lot of instances of scale-free conservative cascades are 
known to evolve steady power-law spectra. Now we can understand how they arise. Indeed, it 
follows from the definition of spectrum (jSJ, the properties (fTTj) and (fT^ . and the arbitrariness 
of S* in (HSl), that 

n{s,so,E) = E^^^X'^n{Xs, Xso,l), for any positive A, 

r = l + — . (16) 
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Hence n has the form: 

n{s,s,,E) = E^"^s-^f{s^/s). (17) 

For wide classes of interactions the right hand side of (|T5|) has a finite hmit as Sq/S 0. 
These classes depend on asymptotic properties of interactions. They have been determined 
for cluster coagulation and weak wave turbulence ('locality conditions') In this 

limit the scale sq drops out of the equations, and the function / has a finite limit: 

lim/(x) = C, and, thus, n(s) ^ CE^/'^s^". (18) 

In turbulence the constant C was first introduced by Kolmogorov 5], and is often denoted as 
6*2. It is clearly interaction-dependent and, therefore, non- universal. For cluster coagulation 



an expression for C was found in 



24|. 



At this point a remark is appropriate: the choice of size variable, which one faces looking 
for a model for a cascade process, is quite arbitrary. For example, one could characterize 
cluster coagulation by the cluster diameter distribution, rather than the cluster mass distri- 
bution. Let the new size variable s = s'^. Then h{s)ds = n{s)ds, and (r — 1) ^ g(r — 1). 
Clearly, (fT^ must obey this transformation rule. It does so, indeed, as under this transform 
{m, a, h} — > {qm, qa, h}. 

The self-similar distribution (|18|) is our central result. Table 1 demonstrates that for many 
complex processes, which can be viewed as scale-free conservative cascades, the exponents 
of power-law spectra can be obtained in nearly automatic fashion. More of examples of 
self-similar spectra can doubtlessly be found scattered throughout the scientific literature in 
various fields. Because values of the indices for a particular cascade are often apparent due 
to their clear physical meaning, our result has considerable heuristic powers. 

The semblance between phenomena in the fields of turbulence, cluster coagulation (aggre- 
gation) and fragmentation, self-organized criticality, critical phenomena, and, more recently, 



econophysics - has long been perceived, e.g. 



30[ . In certain instances it has led to new quan- 



titative results: e.g. the self-similar spectra due to cluster coagulation have been derived for 
three coagulation mechanisms - Brownian, laminar shear, and sedimentation in gravity - 
using reasoning patterned on the Kolmogorov 1941 theory of turbulence Q]. Our approach 
justifies and generalizes such results. 

It is notable that most of well-studied cascade processes - branching processes, fragmen- 
tation, (strange/anomalous) diffusions. Levy flights, multiplicative processes, and cascade 
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models of turbulence intermittency - are linear cascades with h — 1. Although just a par- 
ticular class of scale-free cascades, their investigation has led to rich developments in our 
understanding and modeling of complex phenomena. We expect, therefore, that investiga- 
tion of nonlinear cascades is essential to further progress in this field. 
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